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Abstract 


The Space Shuttle wing-leading edge consists of panels that are made 
of reinforced carbon-carbon. Coating spallation was observed near the 
slip-side region of the panels that experience extreme heating. To 
understand this phenomenon, a root-cause investigation was conducted. 
As part of that investigation, fracture mechanics analyses of the slip-side 
joggle regions of the hot panels were conducted. This paper presents an 
overview of the fracture mechanics analyses. 


1.0 Introduction 

Each Space Shuttle Orbiter wing is comprised of 22 leading edge panels that are made of 
reinforced carbon-carbon (RCC). These panels are part of the thermal protection system that 
protects the wings from extreme heating that takes place during entry into Earth’s atmosphere 
[1], On some of the panels that experience extreme heating, spallation of RCC coating was 
observed in the slip-side regions of the panels. To understand the reasons for this coating 
spallation anomaly, a root-cause investigation was conducted. The root-cause investigation team 
consisted of researchers from various disciplines such as structures, materials, non-destructive 
evaluation, aerothermal, flight hardware, arc-jet testing, coupon testing, vibro-acoustic testing, 
aeroloading testing, photo-micrographic investigation, etc. This paper describes an overview of 
the structural and fracture mechanics analyses that were conducted. 

The paper is organized as follows: First, the Shuttle wing-leading-edge (WLE) configuration is 
presented, and the loading environments to which the Shuttle is subjected are described. Next, 
the building-block approach and analysis models that were used for the stress and fracture 
mechanics analyses are described. Then, the fracture mechanics analyses and methods used to 
characterize defects in the panels are presented. Finally, analysis results and findings are 
discussed. 

2.0 Wing-Leading-Edge Configuration 

The Space Shuttle WLE configuration and geometry are presented in Figure 1. Between 
adjacent panels, T-seals are installed to close or cover the gap. Figure 1(a) shows the Space 
Shuttle WLE, and Figures 1(b) and 1(c) show close-up and cross-sectional views, respectively, 
of Panels 9 and 10 and T-seal 10. The regions of the panels that have S-shaped curvature are 
called “joggles” (see Figure 1(c)). When T-seal 10 is installed, it is locked to a joggle on Panel 
9; that side of the panel is called the “lock side” and is shown in Figure 1(c). In addition, T-seal 
10 floats over a joggle on Panel 10; that side of the panel is called the “slip side” and is also 
shown in Figure 1(c). The region between the lock side and slip side of each panel is called the 
acreage region. The outer mold line (OML) and inner mold line (IML) are the configuration 
boundaries on the outer and inner surfaces of the WLE, respectively. The chord direction runs 
around the WLE and is defined in Figure 1(b), while the span direction runs across the WLE 
panels and T-seals and is defined in Figure 1(c). The RCC material contains a substrate region 
and two conversion coating regions containing craze cracks, as shown in Figure 1(d). Note that 
these craze cracks form in the coating during the cool-down phase of the conversion process of 



fabrication. Then a sealant is applied to both the OML and IML surfaces, filling the craze 
cracks. The conversion-coating regions are commonly referred to as coating layers. 



a) Space Shuttle WLE. b) WLE Panels 9 and 10 and T-seal 10. 



c) Cross-sectional view of Panels 9 and 10 d) Cross-sectional views of slip-side joggle, 
and T-seal 10. 

Figure 1. Space Shuttle WLE configuration and geometry. 

After two different Space Shuttle landings (Missions STS- 102 and STS- 103) at the NASA 
Kennedy Space Center, small areas of the RCC coating were observed to be missing from the 
WLE apex near the top of the joggle in the slip-side region, as shown in Figure 2. The two 
spallation events occurred (a) after space transportation system (STS)- 103 on orbiter vehicle 
(OV)-103 Panel 8L and (b) after STS- 102 on OV-103 Panel 10L. After the Shuttle landed 
following the return-to-flight mission (STS- 1 14), an infrared thermography indication (from 
non-destructive evaluation (NDE)) was observed on OV-103 Panel 8R (Figure 2(c)), indicating a 
loss of coating strength. In addition, the spallation event on the nose cap on OV-105 was 
observed during refurbishment (Figure 2(d)). The root-cause investigation was subsequently 
conducted to determine what factors and scenarios contribute to these spallation events (e.g., see 
References 1-6). 
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Slip Side Joggle Gap 



a) Post STS- 103: OV-103 Panel 8L. b) Post STS- 102: OV-103 Panel 10L. 


c) Post STS-1 14 themography: d) Observed during refurbishment: 

OV-103 Panel 8R. OV- 105 nose cap. 

Figure 2. RCC spallation history. 

3.0 Loading Environments 

During a mission, the Space Shuttle is subjected to four distinct loading environments [1], 

During lift-off and ascent, aerodynamic loads act on the Shuttle; while the Shuttle is in orbit, it is 
subjected to extreme cold temperatures; during entry, the Shuttle experiences peak heating; and 
during descent and landing, different aerodynamic loads act on the Shuttle. 

Loads analyses showed that the lift-off and ascent loading bounds the loading on descent and 
landing. Thus, fracture mechanics analyses were performed using a bounding pressure load over 
an assumed defect for the lift-off and ascent loading condition. The bounding pressure load 
occurs at the maximum dynamic pressure condition (max-Q) about 80 seconds after lift-off. The 
defect driving forces were calculated for the bounding pressure load, and they were nearly zero. 
Hence, both the lift-off and ascent conditions and descent and landing conditions were dismissed 
as being benign contributors to the spallation root cause [6], 

During entry, the Shuttle WLE panels experience entry heating that depends in part on the entry 
trajectory. For a representative entry trajectory from the International Space Station, the peak- 
heating entry temperature along the slip-side joggle regions for each of the 22 WLE panels is 
indicated in Figure 3. Peak heating reaches temperatures up to about 1650°C (3000°F) and 
occurs in the hot panels, Panels 8, 9, and 10. For the on-orbit heating/cooling thermal cycle, the 
temperature range is -130°C < T< 90°C (±200°F), with -130°C (-200°F) being the cold 
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condition. In this paper, fracture mechanics analyses are presented for these two bounding 
thermal loading conditions (i.e., on-orbit cold and entry peak heating) that were performed in 
support of the root-cause investigation. 
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b) RCC WLE slip-side inboard joggle non-catalytic peak temperatures. 

Figure 3. Typical peak entry temperature distribution for the Space Shuttle WLE. 

4.0 Building-Block Approach and Analysis Models 

A building-block analysis approach begins with basic elements and builds in complexity in a 
systematic and progressive manner [2], Such an approach permits each step in the process to be 
verified and its influence on the overall response determined. The hierarchy of the building- 
block approach is presented in Figure 4. First, the problem, objectives, products, resources, 
schedule, and stakeholders are defined. Second, the structural configuration, loads, materials, 
boundary/interface conditions, and tools to be used to solve to problem are identified. Next, 
analysis models, such as finite-element models of the individual components, are created. The 
material modeling procedure, boundary conditions, and other approximations are assigned to the 
individual component models, and the models are solved and verified for accuracy. The process 
at the individual component level is an iterative process that is repeated until the results can be 
verified by comparison to reference solutions or test data. Then, the individual component 
models are assembled, incorporating component interface conditions, different design 
configurations, and various load cases. The assembled models are solved, and the results are 
assessed. The process at the assembly level is also an iterative process that is repeated until 
confidence in the results can be demonstrated and advocated by the analysis team itself. Finally, 
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the results are presented to external reviewers and the stakeholders, engineering reports are 
written, and the entire effort is subjected to rigorous peer review. If necessary, the process may 
be repeated from any of the building-block steps. 



Peerreview 
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Problem definition, objectives, products, resources, schedule, stakeholders 



Figure 4. Building-block approach. 

This building-block analysis approach provides different entry points for analysts, stakeholders, 
and external reviewers. Analysis assumptions can be challenged and revised, analysis findings 
can be compared to observed behavior, and bounding responses can be determined; and the 
analysis end products have technical rigor and their applicability limits defined. This type of 
building-block approach was applied to the global structural and fracture mechanics analysis 
effort as part of the RCC spallation root-cause investigation. All finite element analyses reported 
in this paper were performed using the ABAQUS™ 1 commercial software [7], 

4.1 Analysis Models 

The finite element models used for the analysis effort are presented in Figure 5. First, an 
integrated model was considered (see Figure 5(a)) including Panels 9 and 10 (i.e., two of the hot 
panels), their adjoining T-seal, and attachment hardware. The peak-heating temperature 
distribution was applied to this integrated model. Negligible interaction was observed between 


1 ABAQUS™ is a registered trademark of Dassault Systemes. 
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Panel 10 and T-seal 10 for the peak-heating conditions. These results demonstrate that a single 
panel can be analyzed alone for the entry thermal conditions. 

Next, the Panel 10 model shown in Figure 5(b) was considered because its slip-side joggle region 
experienced the highest entry temperature of the integrated model. The model had 3D elements 
all along the chord in the joggle regions and shell elements in the acreage regions. The peak 
entry temperature distribution shown in Figure 6 was applied to the model. The resulting 
through-the-thickness (TTT) stress distribution in the substrate is presented in Figure 7. (To 
show the TTT stresses clearly, the coating elements are removed in Figure 7.) As shown in 
Figure 7(c), an elevated stress value is observed all along the slip-side joggle in the local region 
near the panel apex (i.e., the hot region) and is nearly constant all along the chord direction (i.e., 
the gradient of stress in the chord direction is nearly zero). These results demonstrate that any 
slice perpendicular to the chord direction could be taken and used in a plane-strain analysis of the 
slip-side joggle [5]. 

Finally, the plane-strain model of the Panel 10 slip-side joggle region (shown in Figure 5(c)) was 
considered. For the stress analysis, the finite element model did not include craze cracks or 
defects. Because the temperatures in the local region near the panel apex are nearly constant (see 
Figure 6(c)), the maximum peak entry temperature was applied to the entire plane-strain model. 
The resulting TTT stress distribution in the substrate (the coating elements have been removed) 
is presented in Figure 8. 

A failure criterion (or criteria) that is (or are) based on linear-elastic stress analysis results alone 
would suggest that wide-spread spallation would occur on many WLE panels on every mission, 
while the flight history shows that wide-spread spallation did not occur. Note that the observed 
spallation events shown in Figure 2 from flight history were limited to a single event at a single 
location on a single panel on each occurrence. However, the stress analysis results are used to 
indicate locations where potential subsurface defects may contribute to a spallation anomaly and 
prompted the following fracture mechanics analyses. 
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Panel 9 


T-seal 10 


Panel 10 



a) Integrated model. 






b) Panel 10 model. 
Joggle Shoulder 



c) Plane-strain model. 
Figure 5. Analysis models. 
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a) Panel 10 model. b) Slip-side joggle region. c) Local region near apex. 

Figure 6. Peak entry temperature distribution for Panel 10. 



a) Panel 10 model. b) Slip-side joggle region. c) Local region near apex. 

Figure 7. TTT stress distribution for peak-heating condition for Panel 10 model. 
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Figure 8. TTT stress distribution for peak-heating condition for plane-strain model. 

5.0 Fracture Mechanics Analyses 

The objective of the fracture mechanics analyses is to evaluate the defect driving forces, which 
are characterized by the strain energy release rates, and determine if defects can become unstable 
for each of the loading conditions. 

The fracture mechanics analyses were performed using the plane-strain model. For the fracture 
mechanics analyses, subsurface defects were introduced at the maximum stress location (see 
Figure 8) and perpendicular to the maximum tensile TTT stress. Three types of defects were 
considered: interface defects (those along the coating-substrate interface), substrate defects 
(those completely within the substrate), and combined interface-substrate defects (those in both 
regions). Typical interface and substrate defects are presented schematically in Figure 9. In this 
paper, only interface and substrate defects are discussed. 



b) Substrate defect. 

Figure 9. Typical interface and substrate defects. 

The finite element model and terminology are presented in Figure 10 (mesh is excluded for 
clarity). The panel acreage is to the left of the model. (Note that the orientation of Figure 10 is 
opposite to Figure 8. The orientation of Figure 10 is used in the remainder of the paper.) The 
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substrate and coating materials are shown in yellow and blue, respectively. The defects are 
inserted by introducing coincident nodes, as demonstrated in Figure 11. To insert a subsurface 
defect between elements A-D and E-H, nodes 2-4 are duplicated to create new nodes 2'-4'. The 
connectivity of elements A-D is left unaltered, while the connectivity of elements E-H is 
modified to use the new nodes 2'-4'. These new nodes are initially coincident with the original 
nodes; i.e., both sets of nodes have identical coordinates. Nodes 1 and 5 define the defect tips. 
Similarly, the craze cracks are modeled explicitly by using coincident nodes in the finite element 
mesh. The craze cracks are oriented normal to the joggle surface, as shown in Figure 10. The 
defects are placed symmetrically and directly underneath a craze crack as specified by the Defect 
Location Index. Location 0 corresponds to the craze crack at the joggle shoulder (see Figure 
5(c)). The craze cracks are numbered consecutively from Location 0. Negative-numbered craze 
cracks are consecutive towards the panel acreage, and positive-numbered craze cracks are 
consecutive down the joggle. The left and right tips of the defect are towards the acreage and 
towards the joggle, respectively. 


OML 


Defect Location 
-2 -1 0 Index 



Figure 10. Plane-strain model for fracture mechanics analyses. 



Figure 11. Defect modeling in the finite element mesh. 

Interface and substrate defects interact with the craze cracks differently for elevated temperature 
conditions. As shown in Figure 9(a), the craze-crack surfaces above an interface defect are 
completely disconnected from the substrate material. As heating occurs, the two craze-crack 
edges may come into contact and try to slide past each other, as shown in Figure 12(a). For this 
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reason, friction is included in the finite element model to account for the craze-crack edge 
interaction. In contrast, as shown in Figure 9(b), the craze-crack surfaces above a substrate 
defect are connected to substrate material. When heating occurs, the two craze-crack edges may 
come into contact, but they are restrained from sliding past each other, as shown in Figure 12(b). 



a) Interface defect with craze-crack edge interaction. 

-1 0 



b) Substrate defect. 

Figure 12. Deflection of interface and substrate defects. 
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The capability of a structure with an embedded defect can be described using the total defect 
driving force Gj calculated using fracture mechanics analyses; i.e., the problem is no longer a 
stress analysis problem. Once the defect driving force is computed, it can be compared to 
material toughness values for the different modes of fracture (Mode I, II, and/or their mixity for 
plane strain). 

In this paper, the G calculation for the defects is performed via the Virtual Crack Closure 
Technique (VCCT). In finite element analysis, the VCCT uses the nodal forces at the node at the 
defect tip and the displacements at the nodes behind the defect tip [8-10]. In the analyses in this 
paper, the finite element models were developed such that the element edge length (A a) of the 
elements ahead and behind the defect tip was of equal size. Figure 13 shows local modeling at 
the defect tip in a 2D analysis. The defect tip is represented by node i. Elements / and J 
contribute to the internal forces, Z, and X h at node i. The internal forces at node i and the relative 
displacements between nodes j and k, A Wjjk = wj - wt and A u^k = uj - Uk, are used to evaluate the 
individual G component values as 

ZiAwjk X i Auj k 

Gt — Git — 

2A a 5 2A a ( 1) 

In Equation 1, G\ and G\\ are the Mode I and Mode II strain energy release rates, the G values. 
The total strain energy release rate Gj can be evaluated using 



Figure 13. VCCT scheme for 4-node (linear) 2D elements. 
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In this paper, a somewhat conservative comparison is made by comparing the total strain energy 
release rate Gj to the Mode I fracture toughness G\ c , which is the smallest of the fracture 
toughness values for the different modes. Unstable defect growth (i.e., defect is likely to grow in 
a catastrophic manner) may occur when Gr is greater than or equal to G[ C ; i.e., GjlG\ c > 1. 

Alternately, a mixed-mode fracture criterion is used to determine whether the defect is stable or 
unstable. In this paper, the 2D plane-strain criterion given by Benzeggagh and Kenane (B-K) 

[1 1] is used as the mixed-mode fracture criterion. This B-K criterion is expressed as 



(3) 


where Gj represents the total defect driving force, given in Equation 2, and R represents the 
resistance to the total defect driving force. The resistance R is also referred to as G c , the critical 
fracture toughness based on test data. The value for the resistance R is given by 


R-G l c + (G)i c - G ]c ) 


( n 

Gn 

\ g t 


- Gjc + (Guc ~ G lc )fi r ' 


(4) 


where Gi c and Gu c represent the Mode I and Mode II material toughness values, respectively, p = 
G\\!Gj represents the mode-mixity ratio, and m represents the exponent of p used to fit the 
mixed-mode test data as shown systematically in Figure 14. When the ratio of Gj to R is less 
than unity, the defect is said to be stable, and defect growth is unlikely. When the ratio equals or 
exceeds unity, the defect is said to be unstable, and catastrophic defect growth is likely. 



g ii/ g t 

Figure 14. Typical curve fit of mixed-mode fracture data. 
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For a mixed-mode response, the B-K resistance R given in Equation 4 is dependent on the value 
of the exponent m as indicated in Figure 14. This exponent is typically determined using a curve 
fit for mixed-mode fracture data obtained from material testing. For a value of m = 1, a linear 
response like the green dashed line in Figure 14 is obtained, which is a non-conservative 
criterion. Mixed-mode fracture tests of polymeric matrix composites show that the B-K criterion 
using m = 2 or 3 appears to fit mixed-mode fracture data accurately. A typical mixed-mode 
fracture criterion curve is shown as red in Figure 14. For mode-mixity ratios less than p, the 
response is dominated by Mode I behavior. 

Single cantilever beam and end-notch flexure test configurations were used in RCC fracture 
testing to determine G\ c and G\\ c . In the absence of further test data, a value of two was selected 
as the value of the exponent m. In addition, the B-basis fracture toughness values are used in 
evaluating Equation 4 as they provide a smaller value for R, which is conservative. 

6.0 Results and Discussion 

In this section, only representative results in terms of defect driving forces, the Gj values, for 
entry peak heating and the on-orbit cold condition are presented. Both interface and substrate 
defects are considered. More detailed analyses and results can be found in References 1-6. 

6.1 Interface Defects - Effect of Defect Location 

Interface defects were studied by assuming a constant-length defect at different craze-crack 
locations. Both entry peak heating and the on-orbit cold condition were considered. The results 
for the normalized Gj values are presented as the ratio Gj!G\ c . Values of Gj!G\ c less than unity 
indicate that the defect is stable and unlikely to grow catastrophically. Values greater than unity 
require further study. Such values could mean the defect is unstable and likely to grow 
catastrophically, or they could mean the mode of fracture is not Mode-I dominated. Hence, the 
total and individual component G values need to be calculated and assessed. 

Entry Peak Heating. For entry peak heating, the craze-crack surfaces come into contact and slide 
past each other, as shown in Figure 15(a), and friction along the craze-crack edges influences the 
defect deformations. Figures 15(b) and 15(c) present the normalized Gj values for two values of 
the coefficient of friction // for various locations of the defect along the joggle. For each 
location, the interface defect is centered beneath the craze crack. Both the left and right tip 
normalized Gj values are plotted as a function of the defect location. For // = 0 (Figure 15(b)), 
the left tip values are considerably smaller than the right tip values for all defect locations 
considered. For both ju= 0 and // = 0.4, the maximum normalized Gj occurs at the shoulder of 
the joggle, Location 0. However, the normalized Gj for // = 0 is considerably higher than the 
normalized Gj for // = 0.4 indicating that the craze-crack edge interaction is an important 
variable to the defect driving force. 

The individual mode G values were also examined. For both ju = 0 and // = 0.4, the G\ values are 
larger than the Gw values. This result suggests that the defect opening for entry peak heating is 
dominated by Mode I. 
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a) Interface defect at craze-crack Location +1 and deformed configurations for interface defects 
at varying craze-crack locations (deformations scaled by a factor of 20). 



Horizontal Distance From Location 0 (craze crack spacings) 


b) Normalized Gj for = 0. 



Horizontal Distance From Location 0 (craze crack spacings) 
c) Normalized Gj for ju = 0.4. 

Figure 15. Normalized Gj as a function of interface defect location - Entry peak heating. 
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On-Orbit Cold. For the on-orbit cold condition, the craze-crack surfaces displace away from 
each other, as shown in Figure 16(a), and the craze-crack edges do not interact and do not 
influence the defect deformations. Figure 16(b) presents the normalized Gj values as a function 
of the interface defect location. The maximum normalized Gj occurs in the acreage area (for the 
left tip at Locations -1 and 0). As the defect is moved into the joggle region, the normalized Gr 
values decrease and reach a plateau. Both the left and right tips have nearly the same normalized 
Gj values except at Location 0, where the left tip (towards the acreage) has the higher value. 

This difference in behavior at Location 0 is due to the geometry; at Location 0, half of the defect 
is in the acreage, and half is down the joggle curve. 

The individual mode G values were also examined. The Gn values are larger than the Gi values, 
suggesting that the defect opening for on-orbit cold is dominated by Mode II and indicating that 
the comparison of Gj to G\ c in Figure 16(b) may be overly conservative. 

-1 0 


-1 0 





a) Interface defect at craze-crack Location +1 and deformed configurations for interface defects 
at varying craze-crack locations (deformations scaled by a factor of 20). 



Horizontal Distance From Location 0 (craze crack spacings) 

b) Normalized Gj. 

Figure 16. Normalized Gj as a function of interface defect location - On-orbit cold. 
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6.2 Substrate Defects - Effect of Defect Length 


The fracture mechanics response for substrate defects is examined by fixing the left tip of the 
defect at different acreage locations and then extending the defect along the joggle (i.e., changing 
the defect size or length). The notation ‘A— >B’ is used in this paper to indicate that a subsurface 
defect exists between Location A and Location B. The left and right defect tips are located at 
Locations A and B, respectively, as illustrated in Figure 17(a). 


Entry Peak Heating. To study “very long” substrate defects for entry peak heating, the left defect 
tip is placed at Location -2, and the location of the right defect tip is varied along the joggle up 
to Location +4. The fracture response is given in Figure 17(b) for the left tip and in Figure 17(c) 
for the right tip. The blue curves represent the total defect driving force, or Gj, normalized by 
Gic, and the green dashed curves represent the resistance to the defect driving force, or R, 
evaluated using Equation 4 and normalized by G\ c . For the initial case (i.e., case —2 — >0), both 
defect tips are stable. For longer defects (i.e., cases -2— >+l and -2— >-+2), the left tip is stable, 
but the right tip is unstable. For the case -2— H-3, the normalized R at the right tip becomes 
larger than the normalized Gt suggesting a return to a stable region. For the “very long” defect 
(i.e., case -2— >+4), the left tip normalized Gt exceeds the normalized R and hence the left tip 
becomes unstable. 
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a) Various defects with the left tip fixed at Location -2; case -2— >+4 shown. 
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b) Left tip. c) Right tip. 

Figure 17. Fracture mechanics response for the left substrate defect tip fixed at Location -2. 


17 


7.0 Effect of Various Variables on the Defect Driving Force 


As mentioned previously, there are many variables that contribute to the defect driving force, the 
Gj value, for this application. Some of these variables are identified in Figure 18, and the effect 
of each of these variables is discussed qualitatively. 


C- 3D effects -JD 



Figure 18. Effect of various variables on defect driving force. 

The effect of the variable in the green cloud (i.e., defect initiation site) was determined by 
studying the photomicrographs like the one shown in the center of this figure. 

The effect of the variables in the orange clouds was determined [1-6] through finite element 
and fracture mechanics analyses like those presented in this paper. For the coating/substrate 
transition zone, the initial analyses assumed a sharp coating/substrate interface and yielded 
higher Gj values. Because a sharp interface does not exist in reality, analyses were conducted 
using different size zones where the material properties were transitioned from substrate to 
coating, yielding lower Gj values. For defect size, larger defects yield higher Gt values. For 3D 
effects, part-through defects yield lower Gj values than those computed using plane-strain 
analysis (i.e., simulating a through defect). These results demonstrate that plane strain is the 
bounding case. For the extent of fiber bridging, including fiber bridging yields lower Gt values. 
For the extent of ply convolutions, including knuckles, voids, etc. yields higher Gt values. For 
stress-free temperature, using a stress-free temperature that increases the AT yields higher Gt 
values. For craze-crack edge interaction, including friction in the analysis affects the results; 
high values for the coefficient of friction yield lower Gt values, while simulating a smooth edge 
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using zero friction yields higher Gj values. For craze-crack orientation, craze-cracks oriented 
normal to the joggle yield the lowest Gr values. 

The effect of the variables in the gray clouds (i.e., defect growth model and nonlinear stress- 
strain response) could not be quantified because enough RCC material was not available to 
perform comprehensive testing, in addition, testing at elevated temperatures is much more 
complex than at room temperature. 

8.0 Concluding Remarks 

The Space Shuttle wing-leading edge consists of panels that are made of reinforced carbon- 
carbon. On some panels that experience extreme heating, spallation of coating was observed in 
the slip-side joggle regions of the panels. Global structural and local fracture mechanics 
analyses were performed on these panels as a part of the root cause investigation of this coating 
spallation anomaly. The global structural analyses showed minimal interaction between adjacent 
panels and T-seals that bridge the gaps between the panels for entry thermal conditions. A 
bounding temperature distribution was applied to a representative panel, and the resulting stress 
distribution was examined. For this loading condition, the through-the-thickness normal stresses 
showed negligible variation in the chord direction and increased values in the vicinity of the slip- 
side joggle shoulder. As such, a representative span-wise slice on the panel was taken, and the 
cross section was analyzed using plane-strain analysis. In the plane-strain models, both interface 
and substrate defects were introduced. Various size defects were considered. Plane-strain finite 
element analyses were conducted for entry peak heating and on-orbit cold conditions. Defect 
driving forces in terms of the strain energy release rates were used to characterize the defects. 
This paper presents some of the fracture mechanics analyses and results. 

Various parameters that affect the driving force of defects present in the Space Shuttle wing- 
leading-edge joggle regions have been identified. In this paper, the effects of defect initiation 
site and defect size were briefly presented. For the fracture mechanics analyses, interface, 
substrate, and combined defects were introduced into the 2D plane-strain finite element models 
of the slip-side joggle region. 3D analyses showed that the plane-strain case is the bounding 
case. The defect driving forces were computed using the Virtual Crack Closure Technique, and 
the Benzeggagh-Kenane mixed-mode fracture criterion was considered. Parameters that affect 
the defect driving force include effects of the coating/substrate transition region, fiber bridging, 
ply convolutions, stress-free temperature, and craze-crack orientation and edge interaction. The 
fracture mechanics analysis results were ultimately used to define tests and test methods and 
helped in understanding possible factors and scenarios that contribute to the RCC spallation 
anomaly. 
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